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Abstract 

In this treatise the well-known 2-body problem with a rotating central body is systematically reinvesti- 
qq 1 gated on the basis of the Projective Unified Field Theory (PUFT) under the following aspects (including 

£f~^ ' the special case of the Newton mechanics): First, equation of motion with abstract additional terms be- 

, ing appropriate for the interpretation of the various effects under discussion: tidal friction effect as well as 

£N| ' non-tidal effects (e.g. rebound effect as temporal variation of the moment of inertia of the rotating body, 

general-relativistic Lense-Thirring effect, new scalaric effects of cosmological origin, being an outcome of the 
. scalarity phenomenon of matter (PUFT). Second, numerical evaluation of the theory. 

Key WORDS: two-body problem with rotating central body - tidal and non-tidal effects - scalaric-cosmological 
■ influence of the expanding cosmos on the 2-body system. 

Q-r 

8 ■ 1 Basis of the new approach 

Caused by the enormous progress in astrophysical precision measuring technique during the last decade (Williams 
et al. 2003, Anderson et al. 2002, Schneider 2003; further references in these papers), mainly based on the 
satellite tracking technique, especially the satellite laser ranging, the theoretical interest in this field of research 
grew considerably. Obviously the main part of the transfer of rotational angular momentum from the rotating 
earth to the orbital angular momentum of the whole earth-moon system has its origin in the tidal braking of the 
mentioned rotation by the gravitational interaction between the orbiting moon and the flowing viscous waters 
of the oceans (Lambeck 1980, Brosche and Schuh 1998, Sabadini and Vermeersen 2004; further references in 
these publications). 

Nowadays particular empirical investigations refer to the additional non-tidal effects, since because of the 
high precision measuring technique, the latter can empirically be separated from the tidal effects. In this 
context some further fields of research are investigated now: Apart from the post-glacial rebound effect already 
mentioned in the above quoted monograph of Lambeck and in greater detail investigated in a recent publication 
(Wu, Schuh and Bibo 2003, further references in this paper), in the general-relativistic Einstein mechanics the 
change of the gravitational field outside a rotating body, particularly the "gravito-magnetic" contribution of the 
Lense-Thirring effect should be mentioned (Lammerzahl and Neugebauer 2001, where further basic references 
can be found). Beyond that the current theoretical and empirical investigations aim to various other scientific 
subjects. Our present field-theoretical research on the basis of PUFT points to new scalaric-cosmological results 
being also interesting for application to the earth-moon system (Schmutzer 2004). 

*eschmu@aol.com 
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1.1 Equation of motion of the 2-body system 

We start with our equation of motion for the case of non-relativistic velocities and weak gravitational and 
scalaric fields (reference mentioned above) and apply it to two different bodies. Then as usual we pass over 
to the combined single equation of motion for the 2-body problem. Having in mind the later application to 
the 2-body system being under cosmological influence, we arrive at the following equation of motion for such a 
mechanical system (dot means derivative with respect to the time t): 

b ) r + ^r + £r + g azlvad = with b) S = 1 ^ , (1) 

where a is the scalaric world function and S is the "logarithmic scalaric world function" , both (of non-Newtonian 
origin) describing the cosmological influence on the 2-body system. Further g az i ta< n is a free abstract azimuthal- 
radial acceleration term introduced for grasping the transfer mechanism of rotational angular momentum from 
the rotating central body to the orbital angular momentum of the 2-body system. 

We remember that (apart from the two additional terms mentioned) the equation Q refers to the frame 
of reference of the center of mass of both bodies. As mentioned above, our treatment of the 2-body system is 
based on an orbiting body 1 (because of thinking of the moon later index M) and an orbiting body 2 (because 
of thinking of the earth later index E), both bodies orbiting about the common center of mass. Of course, the 
theory being developed here can be applied to appropriate 2-body systems occurring in astrophysics in general. 

Let us further mention that 7n is the Newtonian gravitational constant as a true constant of nature. Since 
nowadays there exist serious hints for a time dependence of the empirical gravitational constant G(t), we have 
strictly to distinguish between these two quantities. 

In order to make the reading of this paper easier, we mainly use the notation of our treatment of the 2-body 
system in our textbook (Schmutzer 1989, 2005). 

List of notations: 

Tom (mass of the moon) , (2a) 

rriE (mass of the earth) , (2b) 

M = rriE + tom (total mass of the system), (2c) 

m Em M I A A \ fOA\ 

ni re d — — ~~~ — (reduced mass); (2d) 

ve = — -j^j- r (radius vector from the center of mass to the earth), (3a) 

t*m = — — v (radius vector from the center of mass to the moon) , (3b) 

r = r M - r E , (3c) 
r=\r\. (3d) 

Further following relation is valid: 

m E r E + m M r M = . (4) 

1.2 Conservation of the angular momentum of the 2-body system 

The total angular momentum of the 2-body system reads: 

Ltot — L or i„ + L rot , (5) 

where 

L orb = m red r x f (6) 

is the orbital angular momentum of the 2-body system and 

a) L rot = kL rot with b) L rot = 9 e uj e (7) 

is the rotational angular momentum of the rotating central body {9e moment of inertia, lue angular velocity, k 
constant unit vector orthogonal to the s-y-plane) . 



2 



For physical reasons it seems to be acceptable to demand the conservation law of the total angular momentum 
of the 2-body system: 

*§S-0. (8) 

But one should realize that separate conservation laws do not hold: 

. dL ro t dL rot . dL rot . ■ . . 

a) —^f- = fc dt wltn b ) dt = + , (9) 

M QZ , rad , (10) 



dt 

where M a z\rad is the torque acting on the 2-body system. 

In order to treat the orbital angular momentum of the 2-body system, we first cast the equation of motion 
into the appropriate form 

7n Mm red 

m red r + e r — h m red £r + m red g azlrad = . (11) 

From the equation of motion we find in the usual way the following expression for the torque: 

M az]rad = -m red r x g az]rad . (12) 

Since the torque and the torque causing force F az \ red are related to each other by the formula 

M az \ rad = r x F az \ rad , (13) 

we obtain by comparing this result to the alternative result (|12|) the relationship 

■F az\rad m redOaz\rad ' (14) 

With the help of l(T2")) the balance equation lfTT)|l takes the shape 

J T 

-m red r x g az]rad . (15) 



dt 

1.3 Balance of the mechanical energy of the 2-body system 

In contrast to the demand of conservation of the total angular momentum of the 2-body system, because of the 
irreversible processes (e.g. friction of viscous waters of the oceans of the earth) here we are confronted with a 
non-conservation phenomenon of the mechanical energy. Therefore we have to balance the whole energy content 
of the 2-body system. 

Let us now first define the total mechanical energy of the 2-body system considered: 

+ T rot , (16) 

where the orbital-gravitational energy is defined by 

Eorb\grav — ^orb ~t~ Ug rav (17) 

with 

a) T or t, = — m red r 2 (orbital kinetic energy), b) U grav — — ^ N — — (gravitational energy). (18) 
The rotational kinetic energy reads 

T ra t=-8 E w%- (19) 
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By differentiation from (|16[) results 

dE to t _ dE or b\ grm dT rot . . 

~dT ~ dt + ~dT ■ ( U) 

From these calculations we obtain the following formula for the orbital-gravitational energy of the 2-body 
system, being used in this explicit shape later: 

7? -I , 2 J N M Hired (l> ^ 

-t^orblgrav — ^ Tllred'T ^ ■ l^-U 

Now it is necessary to make some assumptions on the time dependence of several mass quantities. 

In our PUFT (1995) we introduced the concept that the usual mass of a body (being for cosmological reasons 
in principle time-dependent) is the product of two factors: 

m{t) =Ma(t), (22) 

where M. is the time-independent scalmass of this considered body as the primary essential mass (Urmasse) 
being modified by the scalaric world function a of cosmological origin (see Schmutzer 2004) . If this concept is 
accepted, then for m(t) results 

m = mS . (23) 

Application of this concept to the mass quantities introduced above leads to following formulas: 

a) mE = nT>E£) b) rtiM = tiimE , c) M = MS , d) m rec i = m rec iS . (24) 

With respect to the moment of inertia we are confronted with a particular situation: Admitting an additional, 
empirically motivated explicit time dependence of this quantity, we arrive at the relation 

° ° d0 p 

a) E = 9 E S+ Be with b) B = — f . (25) 

at 

For the following it is suitable to split the additional acceleration g az \ ra d occurring in the equation Q into an 
azimuthal and a radial term, using the polar coordinates {i?, ^} in the x-y-plane: 

9az\rad = ^R9rad + e$g az . (26) 

Now using the results l|23|) and l|25|l . after rather lengthy calculations we find the following expression for the 
time derivative of the orbital-gravitational energy: 

dE orb \ grav /l . 2 2j N Mm re d\ , ■ ■ s 

= - y g "W H ) S ~ ™red (R&g a z + RQgrav ) 

= -(T rb - 2U grav )E- m red (R<Pg az + Rg ra d) ■ (27) 
If we apply the abbreviations 

a) dEa ^ rad = m red (RSg az + Rg rad ) , b) = -(T orb - 2U grav )S , (28) 



the expression l|2f)|) reads 

dE orb \ grav dQ s dE az \ rac [ 
dt ~df dt 



(29) 



With respect to the rotating central body, we defined its rotational kinetic energy already above in the usual 
way JTSJ: 

T ro t = \ Oeoje 1 ■ (30) 
Hence with the help of l|25[) for the time derivative results 

dT rot _ dQfnc 1 A 2 _ dQfnc rp y , 1 ° , 2 /n-1 \ 

— -^T + 2 eEUJE ~ dt T ro t£ + 2 Oe^e > (31) 
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where the abbreviation 
dQ f 



I ric 



At 



Oecoe^e < (32) 



for the irreversible heat production of the rotating central body (being braked by the internal friction of the 
viscous fluid component of this body) was used. Let us mention that braking means fulfilling the condition 

LU E < 0. 

Inserting of (|27[) and l|31|) into (|2U|) leads by means of (|28|1 to the final formula 

dE tot _ dQ s dT rot _ dE az \ rad 
dt dt dt dt 



dt 



E u E - {Torb - T rot - 2U gra v) E - m red (R<Pg az + R 9 rad \ (33) 



2 Transition of the equation of motion to cylindrical coordinates 
2.1 Decomposition of the equation of motion 

For practical application of the general theory presented in section ^ it is favorable to use the cylindrical 
coordinates { R, <P, z } . 

Then the equation of motion Q splits into 

R-R<P 2 -\ — ■ ^ .„ + RE + g ra d = (radial equation), (34a) 

(i? 2 + Z 2 ) 3/2 

R@ + 2R& + d>R£ + g az = (azimuthal equation), (34b) 

(R 2 + z 2 Y 

For special use in astrophysics it is appropriate to derive by differentiation of 134b|) the intermediate result 
1 /<Z>\ 2 „ 1 1 . 1 $ 



z + Sz + jj^ = (axial equation). (34c) 



R= 2UJ R -Y^ R + ^ R) -^ RS + SR) -Y^ + 2^ 9 -- (35) 

Inserting into l|34a|) leads to 

R 2 S 2 - 3R 2 - ^ = ~ £ - — {E - E 2 ) + 3RRE + ^±(ER + 3R) - 4 g az + Rg rad . (36) 
it 2<P 2 2^ 2<P 

2.2 General annotation on the adiabaticity approximation 

Since the exact integration of the very complicated system of non-linear differential equations l|34|l is hopeless, 
we use in the following our adiabaticity approximation introduced in our monograph (Schmutzer 2004). 

In astrophysics/cosmology we have to consider the time dependence of physical quantities on two different 
levels of consideration: 

- quick motion of bodies, particles etc. (relatively quick compared to the cosmological temporal changes), 

- cosmological motion (slow time dependence of the cosmological quantities). 

Practically the adiabaticity approximation means that the integrations are performed in such a way that the 
slow cosmological time dependence is neglected, i.e. that during the integration the time plays a non affected 
parameter role. 

2.3 Integration of the axial equation 

This integration is performed under the conditions 

a) 4 < 1 1 b ) R ~ R (37) 
R 
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(near to the x-y-plane and near to circularity). Then from l|34b ) the well-known oscillation equation 

z + Sz + ^-z^O (38) 
results. Using the same abbreviations as in our textbook (Schmutzer 1989, 2005): 



a) C=ir, b) -0=^^, c) u,. = yf^=^-\& t (39) 

and referring to the small friction case (g < wo), then the solution reads (A, B constants of integration) 

z = e -s st [A am(cj z t) + cos(w 2 i)] . (40) 

The exponential factor describes the cosmologically caused damping of the axial motion towards the x-y- 

plane with the damping parameter g. Further uj z means the angular frequency of the orbital motion. The 
corresponding revolution time is 



2tt 



(41) 



whereas for the logarithmic decrement of damping follows 



S. = \3X^- = \t k S. (42) 

Zn+l £ 

2.4 Radial and azimuthal equations of motion in the plane of motion 

In this case both the equations (|34a|l and l|34b|) for z — read: 



a) R-Rn 2 + ?j^- + ER + g rad = or b) g rad = -R + RQ 2 - - RE ; (43) 

a) Rn+(2R+R£)Q + g az = or b) g az = -(Rf2 + 2QR+ RQS) , (44) 

where the abbreviation Q = <P was used. These equations will be exploited later. 

3 Treatment of the conservation of the angular momentum in the 
plane of motion 

3.1 Orbital angular momentum 

In specialization to cylindrical coordinates the orbital angular momentum takes the form 

a) Lorb = kL or i, , where b) £ or 6 = Tn re dF = m re dR 2 Q with c) F = R 2 f2. (45) 

One should remember that in Newton mechanics in the case of vanishing external influence on the orbiting 
system the quantity F (equivalent to the angular momentum) is a constant of motion, but here, as we know 
from l|10ll , because of the rotating central body, L or {, is not a constant of motion. 
Differentiation of (|45b ) leads to the formula 

k^±. (46) 



dt dt 
Comparison to (|10|l gives 

M az \ rad = kM az \ rad = k dL ^ rb , i.e. (47) 
for the corresponding component the relation 

M aArad = ^ (48) 



G 



holds. 

Our next step consists in differentiating of H45b~) and eliminating in the obtained result by means of (|44b ) . 
We find 



° r = -m red Rg az . (49) 
at 

Substituting this expression into (|48|l leads to 

M az \ rad -> M az = -m red Rg az . (50) 

This result coincides with formula (|f 2fl . being obvious if we use the relation l|26[) for elimination. 

From (50) we see that g rad does not enter into the torque, i.e. the torque is determined by g az only. 

Remembering the general relationship between the torque and the torque causing force we arrive at the 
formulas: 

a) F az \ rad -> F az = e$F az with b) F az = ~m red g az . (51) 
Let us remind the reader that we already above calculated the quantity g az (|44b % ). 

3.2 Rotational angular momentum and the conservation law of the total angular 
momentum 

The rotational angular momentum and its time derivative are presented in the formulas J7J and 10. Next we 
write (jHJ by means of and (0i) hi the form 

L tot = kL tot with b) L tot = L orb + L rot . (52) 

Now we insert Q49JI and ^jp) into the conservation law of the total angular momentum (JHJ) in the form 

dLfnf dL nr h dL rn f /r-n\ 

(56) 



dt dt dt 

and arrive at the relation 

m red R 2 Q ftl R \ 1 o 

u E = t. [ — + 2— + S) - w E S - — Oe^e (54) 

UE R ' "e 



or rearranged at 

'J? R^ 



m red R 2 fi + Oe^e 



2 (11 K\ ° 
Oeoje + m re dR ^(^ + 2 ^J + Oeuje 



(55) 



4 Treatment in the plane of motion: energy balancing, temporal 
change of the orbital-gravitational energy, total angular momen- 
tum 

4.1 Balancing of the total energy 

Our calculations in section ^ arrived at the balance equation for the total energy of the system considered 
(|33[) . Further the orbital kinetic energy t\ltth) and the gravitational energy of the 2-body system 1)1 8b) read in 
cylindrical coordinates as follows: 



a) T orb = lm red (R 2 + R 2 tf), b) U grav = _ ^ Mm ™* , ( 56 ) 
Z H 



Hence instead of ()21J) for the plane of motion follows 

T? -T J. 77 - 1 (I?? _L E>2 2\ J^Mm red 

^orb\grav — 1 orb T U grav — ^ m red [tl T-t « ) ~ • (0 1 ) 
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Let us now remember the formulas (|43b ) for g rad , containing R (2nd order derivative of R), and l)44b) for g az , 
containing only first order derivates of the involved physical quantities being interesting for us. Up till now we 
don't have any empirical information on the quantity R. Therefore, of course we have to think about how to 
get some knowledge on this interesting radial second order effect. Here obviously on the way via the quantity 
9rad, i-C by getting further information on this quantity. 

For this reason we cast the balance equation 1|33[) into the form 

dE tot _ dQfnc p v 

— ~^r +EtotS - (58) 

The fulfillment of this postulate is reached by the choice 



Hence the radial equation §43h ) takes the form 



, Q 7nM / R \ 1 

R + RQ 2 + R 2 f2^ + -^- r [l--rE + rfe^| = 0. 60 

R R 2 \ R ' m red R 

Future precision experiments on radial second order effects could show whether this equation is applicable. 
Further we derive with the help of (|54|l from 144b ) 

9az = rf E (Ve + u E £ + -!- Oe^e) (61) 
and from l|28b) the result 

aE az\rad ( t>1 r& , p2 , 1^ M \ . 1 ° 2 fRO\ 

— = -m red \RS2 +R H —)E+-9eu e - ( 62 ) 

Then 1|29[1 goes over to the simpler formula 

a) = E orblgrav S - - Q e lj% or b) — ^ = 9e^ 2 e ■ (63) 

at la 2 at \ m red J 2m red 

4.2 Temporal change of the orbital-gravitational energy 

Then integration of l|63b ) leads to the following expression for the time-dependent orbital-gravitational energy: 

o ^ 

a) E orb i„ rav — — m red A H — Cm red with b) A = - j E E dt (C constant of integration). (64) 

2 2 J m red 

We explicitly learn from these relations that primarily the origin of the increase of the orbital-gravitational 
energy of the 2-body system considered is caused by the time-dependent decrease of the moment of inertia of 
the rotating body (rebound effect). 

Next by means of (|64h ) we derive from (|57J) the equation 



a) tf + R2tf- 2 J^ = C-2A 1 i.e. b) R = + Jc + ^ - WW - 2A . (65) 
R [ ' \ R 

4.3 Explicit conservation law of the total angular momentum 

In order to treat this task we give the equation l|54|) by means of (|25b ) the shape 

a) m red ^-(R 2 Q a ) + E ^-(w E a )+ 9e^e v = or b) (m red R 2 f2 + 9 e ^e) = . (66) 
dt v ' dt v ' dt v ' 

Integration leads to the conservation law of the total angular momentum 

Jo = Ltot = Lorb + L rot — m red R 2 f2 + 9e^e (67) 
(Jo constant of integration). 
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5 Transition to the Keplerian conic motion by using the adiabaticity 
approximation 

Eliminating in (|65b ) the angular velocity fl derived from (|67|l in the form 

n Jo - Oe^e ,„ q , 
n = m~ > ( 68 ) 

we arrive at the equation 



fl2+ (J -W_2^M =g _ 2A 
mi ed R 2 R 



to be treated further now. Since in the parameters M, m re d, Oe, l^e, A being involved in this non- linear 
differential equation a weak adiabatic (partly adiabatic-cosmological) time dependence occurs, we approximately 
take them as adiabatic constants, when we apply the transformation 

R ® = W) (70) 

in order to reach our goal of the adiabatic Keplerian conic equation of motion (e excentricity) 

a) d^ + ^T = Wlth b) T = (j -e EUE y (71) 

The solution of H71fa ) reads in an appropriate form 

R = = . (72) 

A (1 - ecos^) 

In the following we are only interested in the elliptic motion which according to the notion in our monograph 
(Schmutzer 1989, 2005) is determined by (a, b semi-axes) 



e = Va 2 - b 2 < 1 (b < a) . (73) 
Calculating the energy of the 2-body system by inserting of l|68() . 17211 etc. into l|57() leads to the result 

_ 7NMm red 

^orblgrav \ ) 

which formally reminds us of the situation in the Keplerian 1-body system, but one should not forget the 
adiabatic time dependence in the physical quantities occurring in (74). 

Next we calculate from the adiabatically valid formula for the area of the ellipse 

A M = — (75) 

A 2 {l-e 2 ) 3 / 2 

and from the also adiabatically valid Keplerian area law 

dA e ii _ A e ii _ J - Be^e 
dt T e ii 2m red 

the revolution time r e u of the elliptic motion: 

r e u = ~P= a 2/3 • (77) 
Using the angular revolution frequency 

Veil = — , (78) 

Tell 



(76) 
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we find the differential relationship 

dfleii 3 da 
f2 e u 2a 



(79) 



which because of (|74|l is valid for the change of Sl e ii, understandable through the change of the orbital- 
gravitational energy E or b\ grav , caused by external energetic influences. In this context we remember that 
without such external influences (apart from the weak cosmological effects mentioned) and without a temporal 
change of 9e the orbital-gravitational energy according to (|53T would be conserved. 
Let us here remind that the relation l|79|l should not be confused with the relation 

which follows from the conservation law of the angular moment (RrQ = const) of the non-disturbed Keplerian 
1-body system. See equation (45b) for the case F = const. 



6 Influence of the Lense-Thirring effect on the 2-body system 

It is well-known that the Newtonian gravitational field equation is at most usable in the quasi-static case of time 
dependence of the source. Already in the case of a stationarily rotating source the Einsteinian gravitational 
theory leads to a modification of the Newtonian external field which was calculated in first approximation by 
J. Lense and H. Thirring (1918). Nearly half a century later R. P. Kerr succeeded in finding an exact exterior 
solution of a rotating body. 

It is obvious that the Lense-Thirring modification of the gravitational field soon became interesting also 
for astrophysicists. Parallel to the field modification also the modification of the corresponding equation of 
motion of a body was derived by these authors, leading to the Lense-Thirring additional terms in the extended 
Newtonian equation of motion 

* + e ^ " j£ * * " ^ E e r )e r ] = (81) 

(Le — kOEtOE angular momentum of the rotating central body, c velocity of light). A new derivation of this 
equation was recently published (Lammerzahl and Neugebauer 2001, further references in this paper). 
In analogy to the Coriolis acceleration 

a c = 2f ' xw (82) 

one is inclined to define 

w LTh = -^-j [L E - 3(L E e r )e r ] (83) 

as the Lense-Thirring frame dragging angular velocity. 

With the aim to compare the equation l|81|l with our ansatz equation we on the basis of the Newton 
mechanics {S — 0) specialize the equations l|81() and to the plane of motion (e^ and e$ usual unit vectors 
in this plane) and find by using of 126|) : 

R+eR E?{ M ~2 j + e<g C 2 R 3 =0 ' ( 84& ) 

R + e R r2 + e R g rad + e,pg az = . (84b) 
Comparing these two results, we arrive at the identifications 

27N#EW_Ei? _ 2-f N E U! E R 
aj 3rad|LTli , 0) 5 QZ |LTh - ■ (80 J 

In this special case (E = 0) from (|59|l and 144b ') follows: 

a) g Ta d= (2+^l)Rn 2 + —^— E u%, b) g az = -{Rf2 + 2QR) . (86) 
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Of course, the pairs of both the equations H85|) and (|86[) have a fully different physical basis. 

Numerical evaluation shows that in the case of the earth-moon system the order of magnitude of the general- 
relativistic Lense-Thirring effects (18511 are roughly 14 orders smaller than the above calculated effects. Nev- 
ertheless it may be interesting to look how for appropriate star models (yet to look for) the situation could 
change. 

7 Survey of the exterior spherical harmonics expansion of the New- 
tonian gravitational potential including various applications up to 
the second order 

7.1 General mass distribution 

On the basis of the Newtonian gravitational field equation 

A0n = 47T7NAi (87) 

(</>n Newton potential, /i mass density of the source) the expansion considered reads (Lambeck 1980, Schneider 
1999, Schmutzer 1989,2005): 



a < q ^ ^ nM 



r 



1 + E E ( ~ ) P™ (cos z9) <^ C nm cos mip + S nm sm(rmp) > 



n— 1 m— 

where following notation is used: 

M mass of the source body, 

r e radial parameter, 

ip geographic length, 

i) polar angle (0 = — + $ geographic latitude), (89) 

Cn m , S nm Stokes coefficients (tesseral gravitational moments, geopotential harmonic 
coefficients), 

P™ associated Legendre polynomials. 

Choosing the center of mass as origin of the spherical coordinate system, one finds the zero values: 

C w = 0, C u =0, S u = 0. (90) 
Then the equation reads 



9N{r,v,(p) = 



1_ E [-) J « p «( cos ^) 



n— 2 m—1 

where the abbreviation 



EE - P ?T( C0S, $)\ Cnm cos(m(p) + S nm sin(mip) V 



(91) 



J n = —C'no (zonal gravitational moments) (92) 
is used. 

7.2 General mass distribution in second order approximation 

In this approximation, presupposed for our further investigations, the dynamic form factor (quadrupole moment 
parameter) 

•h = -C 20 (93) 
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is of our specific interest. 

For some further calculations it is convenient to go over to Cartesian coordinates (Greek indices run from 1 
to 3). Then the Newton potential takes the form 



^(r)=4 M) (r)+4 w (r) 



with 



TnM 



^( r ) = -^E{ZM aP x a xp - M aa r 2 ) 



where 



M = J n(f) d V 

v 

M af3 = M (ja = I ii{f)x a Xj3 d V 



(monopole potential), 
(vanishing dipole potential), 
(quadrupole potential) , 

(mass of the source), 
(mass tensor of the source). 



By means of this quantity the mass quadrupole moment tensor 
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a) m aj 3 = mpa = 3M a p - M 77 <5 Q/3 with b) m aa = , c) M a p = - (m a p + M^8 a p) 

is defined. Then the potential (O with the help of l(55k). l(9"5bl and (EH reads 
7nM 7 N 



^n(t-) 



2r 5 



m a (3X a Xf3 . 



(94) 

(95a) 

(95b) 
(95c) 

(96a) 
(96b) 

(97) 
(98) 



7.3 Inertia moment tensor and mass quadrupole moment tensor 

The inertia moment tensor 



'a/3 = 6p a = / K r ) (xjXjSa/3 - x a xp) d V 



(99) 



which is related to the mass tensor (|96b ~l and to the mass quadrupole moment tensor H97fa ) via the relations 
a) 9 a p = M 77 (5 Q/3 - M a p with b) 6» 77 = 2Af 77 (100) 

and 



9 Q /3 — -^{p il 5 a p — m a p) 



(101) 



In main-axes representation of the inertia moment tensor the following relationships with the Stokes coefficients 
hold (Schneider 1999): 



a) C 2 o 
a) S 21 



J 2 



'33 — 2 V Pl1 ' ° 22 / 



Mrl 



Mri 



0. 



b) C21 
b) S22 - 



C13 

' Mr 2 e ' 

612 
~2Mrl 



0. 



C) C 2 



AMrl 



Introducing for simplicity the usual notation 

a) 9 n =9 x = A, b) 9 22 = 9 y = B, c) 



733 



= Z =C 



(102) 
(103) 

(104) 
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(A, B equatorial moments of inertia, C polar moment of inertia), we arrive at the following shape of both the 
non- vanishing coefficients: 

a, *-^t£>. -» 

Further notations used in literature are: 

9 aa = A + B + C (trace of the inertia moment tensor), (106a) 

C-UA + B) . . . 

ri = — (dynamic oblateness). (106b) 

7.4 Rotation-symmetric specialization 

Rotation symmetry with respect to the z-axis means for the Newton potential the condition 

£-0. (107) 

Hence for the Stokes coefficients results 

C nm , S nm = for n = 2 , 3 , . . . , (108) 
i.e. from 191fl we find 

7N M 



oo / \ n 
n=2 > r ' 



(109) 



(r e — ► ro equatorial radius of the rotation-symmetric body) . 

Passing over to the second order approximation, we get from l|109f) the expression 

Mr, *) = + (1 + 3 cos(20)) . (110) 

With respect to the mass distribution the rotation symmetry means the conditions 

a) m xx =m yy , b) m xy = m yz = m zx = , i.e. c) m zz = -2m xx = -2m yy . (Ill) 
Now we define the mass quadrupole moment, taking following sign convention: 

m (Q) = -m zz . (112) 
Then, using spherical polar coordinates the formula (98) goes over to 

fo(r,fl) = + 7N 8 T " 3 (Q) (1 + 3cos(2^)) with (113a) 

r 2 =i? 2 + z 2 , i? = rsin^, z = rcost9. (113b) 

Here we would like to mention that the definition of some basic notions in this field of research is rather confusing 
(choice of different prefactors and signs), caused by both the two different historical approaches: nuclear physics 
(see Landau/Lifschitz series of textbooks) and geophysical and general-relativistic research (see Misner, Thorne 
and Wheeler 1973). 

In the rotation-symmetric case applied here the following formulas hold: 

a) e lx = e Xi = ~m^ + ^M ri) b) 033 = |m (O) + §M Tr , c) ia = a3 = 6 3 i = . (114) 

We mention that in this context a basic role plays the polar moment of inertia 

C = 33 = [»(r)R 2 dV (R 2 = x 2 +f), (115) 
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Furthermore, we notice that in the notation l|l(J4l) the formulas 1)1 14b ). (1114b ). (|105fl and 1)106)1 read: 
a) A = B = -±wW + \m t1i b) C =\m^ + \m^ , c) J 2 = , 

d) C 22 =0, e) aQ = 2A+C, f) W f 5 



(116) 



C 



whereas the comparison of 1)113(1 with 1)110(1 gives 
_ m W! _ C* - A 

j2 -2Mnr^p (117) 

Further results are: 

a) M 7J = A+ — , b) m (Q) = 2(C-,4). (118) 

Applying the above outcome to a rotation-symmetric ellipsoidal body with the semi-axes a, c, (a > c) we define 
the geometrical oblateness as usual: 

A ^ = ^ (119) 

Let us now consider this body stationarily rotating about the symmetry axis, where the interior centrifugal 
potential is given by 

(t>cf{r,$) = -\u 2 R 2 = -iwVsin 2 t? = \ u 2 r 2 [l -P 2 (costf)] (120) 
(uj angular velocity, P 2 Legendre polynomial). 

7T 7T 1 

At the equator (e), where r = a, 1? = — , P 2 (cos — ) = P 2 (0) = — - is valid, results 

c /| e = -^ 2 a 2 . (121) 

On the other hand the Newton potential ((110)1 reads on the equator 

7 N M tnM 

•pNle = ; 7^— >h ■ (122) 

a 2a 

Superposition of both the potentials gives the constant total potential on the equator 

7nM / 1 uj 2 a 3 \ 

4>tot\e = <?N|e + <Pcf\e = ~ i 1 + 2 ^ + 2j N M ) = C(MSt ' 

since because of stability the total potential on the surface of the ellipsoid has to be constant. 

Analogous considerations with respect to the north pole (p), where r = c, d = 0, P 2 (cos0) = P 2 (l) = 1 
holds, lead to 



D tot\p 



7nA/ u_(&\ 2 j t 

c L \cJ 



(124) 



Equating the results ()123H and ()124)l . we find 



J 2 = 

'vfl 



^ C L0 2 C 3 

2 a 27 N A/ 



ay v 2 r, M . (i 25 ) 

Using the approximation 

a) A «l, b) (126) 
we arrive at the interesting relationship (see Sharkow 1976) 

between the geometrical oblateness and the dynamic form factor J 2 . 
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8 Tidal and non-tidal influence on the rotation of the central body 

8.1 Splitting of the angular velocity and angular acceleration of the rotating body 

The angular velocity of the central body lje may split into the constant undisturbed angular velocity u>E\ot the 
tidal angular velocity u>E\t an d the non-tidal angular velocity u> E \nt'- 



Temporal differentiation gives the angular acceleration in the form 

LOE = <l>E\t + ^E\nt j 



(128) 



(129) 



where u> E \t 1S the negative tidal acceleration (tidal deceleration because of the braking effect of the oceanic 
waters) and u>E\ n t the non-tidal acceleration which usually is identified with the rebound effect mentioned 
above (positive sign). Our concept has to be enlarged, since according to our aim we also have to grasp the 
cosmological influence. Therefore we separate as follows: 



^E\nt — ^E\expl + ^E\cos 



(130) 



The first term u>E\expi refers to the explicit time dependence of the polar moment of inertia of the central body, 
including the above mentioned rebound time dependence, but admitting further time-dependent effects. We 
call this amount "explicit temporal effect" which according to this interpretation splits as follows: 



^E\expl — ^E\rb + ^>E\further 

The second term d>E\further keeps open the investigation of further effects. 

8.2 Treatment of the explicit temporal effect 

First we start our investigations by temporal differentiation of (|117fl : 



(131) 



Mr* 



C-A 



'"0 



(C - A)r Q 



(132) 



As shown by previous research, in good approximation the earth behaves like an incompressible body, being 
interpreted as a material with constant trace of the inertia moment tensor: 9 = 8 aa = const. Application of 
this idea to l|116b ) leads to the equation 



Inserting into l|132|l gives 
3 



J 2 = 



2Mr 2 



C-2 C- 



( C -H 



If the inequality 

\ 3 / r 



holds, then from (|134f) follows 
3C* 



■h 



2Mrl 



Eliminating C by means of the relation (|25b ) and giving the resulting relation the form 

o o dC 

a) C=C+CS with b) 0= — , 

at 



(133) 



(134) 



(135) 



(136) 



(137) 
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leads to 



j 2 = JLs{i> + CE). (138) 

Let us remember that the tidal braking effect of the earth is an extreme long-term process (billions of years), 
whereas the post-glacial rebound phenomenon is, cosmologically considered, a short-term process (some ten 
thousands of years). Therefore one is tempted to treat (under the relatively constant long-term tidal background) 
the rotating earth as a rather isolated system with an approximately conservative angular momentum (similar 
to the well-known pirouette effect): 

' {Cue) =Cu e +Cw e = 0, (139) 



i.e. 



Of 



° C we 

C = . (140) 

UlE 



Inserting this expression into H138|) leads to the formula 



= -E). (141) 

If one neglects the cosmological influence (£) and reduces the sight only on the rebound effect, i.e. ll>e = w E \ r b, 
then from l|141|) follows 

SCiO E\rh 

J2\rb = - », r % — ■ 142 

2{rb 2Mrluj E 

This formula, on a different basis derived, was recently numerically with success tested for the non-tidal rebound 
effect of the earth (Wu, Schuh and Bibo 2003). 

Apart from the rebound effect, among the further effects being under consideration the hypothetical expan- 
sion of the earth will be shortly treated now. In contrast to the rebound effect with its temporally decreasing 
moment of inertia (increasing angular velocity) , obviously an expanding earth would exhibit an increasing mo- 
ment of inertia (decreasing angular velocity similar to the braking behavior of the earth). For simplicity we 
consider a rotating homogeneous sphere with the moment of inertia 

C = \ Mr'l . (143) 

Considering the expansion process (ep), by explicit temporal differentiation we receive 



2C(dro\ 



C e p = — ( ) >0. (144) 
r V at Jep 



Inserting into 1)138(1 gives the relation 

J2|ep " Mrl 



\( — ) H 




LV dt Jep 





(145) 



Some years ago we hypothetically developed an expansion theory of celestial bodies on the basis of our PUFT 
(Schmutzer 2000, Schmutzer 2004) which led us to the following result for the radial expansion velocity: 

f dr 
\ dt 



= AE s r a 4 T, with (146a) 

ep 

= 7 - tp I (scalaric cosmic expansion factor), (146b) 

f qOLq 

(scalaric material factor), (146c) 



3cc 
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where following quantities are involved: Fq — a \r x f| modified areal velocity of the orbiting body (here moon), 
a c cubic thermal coefficient of dilation of the earth, cq specific heat of the earth, fs heat consumption factor 
of the earth. Let us emphasize that the model for the earth (and therefore also the values occurring) is very 
rough. Nevertheless we pointed as an example to such a further possibility of explicit temporal effects. 

Inserting l|146b ) into (|145|l leads us to the final formula for the temporal change of the dynamic form factor: 

12CE s a 4 T, ( 1 \ , , 

J ^ = ^^( 1+ 8^) >0 - < 147 > 

9 Specialization of the above developed theory to the case of the 
Newtonian mechanics 

In order to simplify our above results based on PUFT, for usual application on the basis of Newtonian physics 
we have to refrain from the cosmological effects (acting on the 2-body system) by setting 

S = 0. (148) 

Then the store of the basic specialized formulas (getting from those listed above) reads: 



a) r + r + e R g rad + e$ g az = , b) R - Ril 2 + + g rad = , 

r 6 R z 

■ 2m i 

c) n + — +K9^ = °> 



^tot — -&orb\grav ^rot j t>) -&orb\grav ~~ ^orb ~r~ Ug rav , c) T or ^ — Tn rec [{^R "f^ ^ ) 
A\ TT jNMmred , rr 1 a 2 



a) 



R ' ' ' 2 

dE ro f ^-&orb\grav dT ro 

dt ~ dl dt 



(149) 



9az = —^—(uj E + — Oe^e) =-(RQ+ 2ttR) , (150a) 



Rrrired v 6e 

9rad = ( 2 + ^2) R^ + - J— E J% , (150b) 

V QRJ 2m red 

m red R 2 n (tl R\ 1 o 

ue = ^— (jj + 2- j - - few* = , (150c) 

R + RQ 2 + ^ + i? 2 fi^ + — ^ = ; (150d) 

a) ifot = L or b + L rot , b) L or b = m red R 2 Q , c) L rot — Oe^e ] (151) 

a) — ^- = , b) — — = M az = m red R^tt^- + 2—j ; c) — j— = 9 E u E ] (152) 

a) M az = F az R = -m red Rg az , b) F az = -m red g az ; (153) 



(154) 



b) — Q ^l9 r< ™ =~\ Qeu e > i- e - E orb\ g rav = -m red A + i Cm^d , (155) 

dT rot dQf r i c 1 ° 2 j\ dQ/ r i c . di? tot dQ/ r ic 

c) ^T = ^7-+2^- d ) - S - = ^W<0, e) — = — <0; 

a) E orb \ grav = „ 7N ^™ red ; L tot = J = m red R 2 fl + 8 E ^E; (156) 

1 dfi e n 3 da 
f2 e ;; dt 2a dt 



(157) 
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10 Numerical evaluation of the theory 

10.1 Main empirical values of physical and astrophysical quantities needed (adapted 
to the 2-body system) 

7n = 6.673 • 10~ 8 g _1 cm 3 s~ 2 (Newtonian gravitational constant). (158) 

Earth (for new values see Morrison and Stephenson 2001, Sabadini and Vermeersen 2004): 

r E = 6.378 • 10 s cm (equatorial radius), (159a) 

m E = 5.976 • 10 27 g (mass), (159b) 

lu e = 7.292 ■ 10~ 5 s _1 (angular velocity), (159c) 

Cjeh = —6.15 ■ 10~ 22 s~ 2 (tidal angular acceleration), (159d) 

&E\nt = 1-35 ■ 10~ 22 s~ 2 (non-tidal angular acceleration), (159e) 

u>E\emp — &E\t + ^E\nt = —4.8 • 10~ 22 s~ 2 (empirical angular acceleration), (159f) 

C E = 8.0394 • 10 44 gem 2 (polar moment of inertia), (159g) 

Ae = 8.0131 ■ 10 44 g cm 2 (equatorial moment of inertia), (159h) 

Ce — —1.2675 • 10 27 gcm 2 s _1 (derivative of the polar moment of inertia), (159i) 

J 2 = 1.0826 • 10~ 3 (dynamic form factor), (159j) 

J2\rb = —1.014 • 10 _18 s _1 (derivative of the dynamic form factor). (159k) 

Moon (for new values of radial velocity and orbital angular velocity see Anderson et al. 2002, Williams et al. 
2003): 

m M = 7.348 • 10 25 g (mass), (160a) 

VM\orb — 1-023- 10 5 cms _1 (orbital velocity), (160b) 

v M\rad = 1-201 ■ 10~ 7 cms -1 (radial velocity) , (160c) 

il M = -1.251 • 10~ 23 s~ 2 (orbital angular velocity). (160d) 

Earth-moon system: 

R = 3.844 • 10 10 cm (distance of the centers of mass), (161a) 

M = m E + m M = 6.047 • 10 27 g (total mass), (161b) 

m re d = mE J jJ JlM ~ 7-259 ■ 10 25 g (reduced mass). (161c) 

10.2 Convenient parameters for fitting the theory to the empirical values 

For simplifying this numerical fitting procedure we introduce 4 adapting parameters a , [3 , 7 , S being under- 
standable from the following equations: 

a) Z' = a-10- 20 s- 1 , b) C E = P ■ 10 44 gem 2 , 

(162) 

c) tb E \emp = -7 • 10 22 s 2 , d) Ce = -5- 10 27 gcm 2 s 1 . 

Since here is not printing space enough for a detailed presentation of the calculations, we only may report that 
the parameters 

a) a = 10.612, b) /3 = 8.0394, c) 7 = 4.8, d) (5 = 1.2675 (163) 

in context with the formulas (|162|) fulfill the decisive equations listed above. The last three parameters directly 
coincide with the empirical results i|159b ). (|159t "). (11591 ). 
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The first parameter follows from the relation 1)550. Inserting this value into H162b ~) leads to the present numer- 
ical value of the empirically important cosmological quantity "logarithmic scalaric world function" introduced 
above, which describes the PUFT-caused cosmological influence on the 2-body system considered: 

£ = 1.06 ■ 1CT 19 s- 1 = 3.35 ■ 1CT 12 y" 1 . (164) 

This result is for our PUFT insofar of basic nature, since this quantity determines the cosmological expansion 
influence on many universal cosmogonical and astrophysical processes. Recently by applying our theory to a 
closed homogeneous and isotropic cosmology without the cosmological term (using the empirical WMAP- values 
of 2003) we arrived at the result 

S = 5.2 ■ 10~ 19 s- 1 = 1.64 ■ 10~ n y" 1 . (165) 

The new value l|164|) is nearly one order of magnitude smaller than the previous value (|165fl . This reduction 
improves the situation with respect to the empirical results of the temporal change of the empirical gravitational 
constant (Schmutzer 2004). 

Concluding this section let us make two further annotations: 

With respect to appropriate future accuracy measurements in this field of research there is good hope for 
an improvement of the measured values by some orders of magnitude in the next years (Williams et al. 2003, 
Turyshev et al. 2003). 

In section El we specialized our PUFT results to Newton mechanics. If one repeats the calculations on this 
basis with S = 0, then one is confronted with some small discrepancies which may be caused by some lack in 
the accuracy of measurements (i.e. Newton mechanics may be correct), or by the mentioned improvement of 
the accuracy these discrepancies are true (i.e. a new theory as e.g. PUFT has to be tested). 

10.3 Report on further numerical results 

Various angular momenta: 



a) L orb = 2.857- 10 41 gem 2 s" 1 , b) L rot = 5.862 • 10 40 g cm 2 s -1 , 

c) L tot = J = 3.441 • 10 41 g cm 2 s" 1 ; 



(166) 



Further: 



a) g rad = -2.609- 10 _1 cms" 2 , b) g az = -1.692 • 10~ 13 cms" 2 ; (167) 
a) 5 rad|LTh = -1-565 -10- 14 cms" 2 , b) g az \ LTh = 1.838 • HT 26 cm s~ 2 ; (168) 

a) F az = 1.228- 10 13 gem 2 s -2 , b) M az = 4.721 • 10 23 gem 2 s" 2 ; (169) 

fgrav 

o ooo m35 „ _2 ..-2 . 

-^orblgrav 

a) T rot = 2.137- 10 36 gem 2 s~ 2 , b) E tot = 1.755 • 10 36 g cm 2 s~ 2 ; (171) 



a) T orb = 3.798 • 10 35 gem 2 s~ 2 , b) U grav = -7.62 ■ 10 35 gems -2 

c) E orb \ grav = -3.822 • 10 35 g cm 2 s~ 2 ; 



(170) 



a) E orb \ grav = 3.329 • 10 18 gem 2 s -3 , b) f rot = -3.128 ■ 10 19 g cm 2 s~ 3 , 

c) Q Sr i C = -2.814- 10 19 gem 2 s" 3 , d) E tot = -2.795 • 10 19 g cm 2 s~ 3 w Q fr 



(172) 



\ d-Qcircle , 3Q c ircle « . n \ ^^circle , n ^circle , n. / 1 r 70 \ 

a) ~AT> + TR ^ b) Tr + 2 R (l73) 

U- r *'circle ^^circle ^-^circle ^circle 

With respect to both these equations we remember that the first one coming from (|79|l is well fulfilled (to be seen 
in connection with the conservation of the total angular momentum), whereas the second similar one resulting 
from H80JI is not fulfilled because of the non-conservation of the orbital angular momentum of the 2-body system 
(in contrast to a situation without external influence). 

I would like to thank very much for scientific discussions P. Brosche, M. Schneider, H. Schuh, U. Walzer and 
Bu Win, and further particularly A. Gorbatsiewich for scientific discussions and technical help. 
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